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ABSTRACT 

We report the first results of an observational program designed to determine the 
luminosity density of high redshift quasars (z>5 quasars) using deep multi-colour 
CCD data. We report the discovery and spectra of 3 i<21.5 high redshift (z>4.4) 
quasars, including one with z>5. At z=5.17, this is the fourth highest redshift quasar 
currently published. Using these preliminary results we derive an estimate of the 
Mb<— 25.0 (Mabi450<— 24.5) quasar space density in the redshift range 4.8<z<5.8 of 
3.6±2.5x 10~* Mpc~^. When completed the survey will provide a firm constraint on 
the contribution to the ionizing UV background in the redshift range 4.5—5.5 from 
quasars by determining the faint end slope of the quasar luminosity function. The 
survey uses imaging data taken with the 2.5m Isaac Newton Telescope as part of the 
Public Isaac Newton Group Wide Field Survey (WFS). This initial sample of objects 
is taken from two fields of effective area ~12.5deg^ from the final '--^lOOdeg^. 
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1 INTRODUCTION 



Tlie identification and study of quasars at liigh redshift pro- 
vides key diagnostic information on the early Universe. Over 
the last few years the relationship between high redshift 
quasars and local galaxies has been given an impetus from 
the realization that quiescent black holes may lie at the cen- 
tre of all local bulges including the Milky Way (e.g. Magor- 
rian etal, 1998). The formation and evolution of massive 
black holes is fast becoming a main stream topic in theories 
of gal axy fo rmation and evolution (Haeh nelt, Natarajan and 
Rees ( |1998| ); Kauffmann and Haehnelt ( ^000| )). 

High redshift quasars are also an important contributor 
to the UV radiation field in the early Universe. Specifically 
the lack of an optically thick, Gunn-Peterson trough in the 
spectra of high r edshift objects (e.g. Songaila, Hu, Cowie, 
McMahon (1999)). implies that th e Uni verse is highly ion- 
ized at z~5.5. Estimates by Madau( 1999) indicate that if the 
ionization is due to starlight the inferred star formation rate 
at z=5 must be comparable to or greater than the observed 
value at z=3. The direct detection of galaxies at z~3 via the 
Lyman Break technique and the detection of measurable Ly- 
man continuum photons that e scape these galaxies (Steidel, 
Pettini and Adelberger (2001)) has raised a question over 



trough should be present by z~6. A major uncertainty in 
all of this work is a measure of the ionizing background 
from quasars. The current uncertainty in the contribution 
of quasars to the UV background is predominantly caused 
by the unconstrained faint end of the luminosity function. 

The resent discovery of several z> 5 quas ars by the Sloan 
Digital Sky Survey (SDSS Fan et al. ( |200l| )) will constrain 
the bright end (i<20.0) of the high redshift luminosity func- 
tion. We are undertaking a survey 1-2 magnitudes fainter 
than the SDSS program with the aim of constraining the 
space density and evolution of lower luminosity quasars 
where the majority of the ionizing flux arises. In this paper 
we present results from the first ~10deg^ of our survey. Un- 
less stated otherwise we use conventional Vega magnitudes 
and Ho =50 km s"^ Mpc"\ qo=0.5. 



2 SURVEY IMAGING DATA 

We are using data from the Isaac Newton Telescope (INT ) 
Wide Angle Survey (WAS) (McM ahon et al. ( ^000| ); 
http://www.ast.cam.ac.uk/~wfcsur/. The survey is car- 



ried out with the prime foc us W ide Field camera (WFC; 
Ives,Tulloch and Churchill (|l996[)) at the 2.5m INT. The 



whether quasars or galaxies are the origin of the UV pho- 
tons that re-ionized the Universe at z>6. Re cent t heoretical 
work (Miralda-Esc ude, H aehnelt and Rees (1998); Madau, 
Haardt and Rees (1999)) indicates that a Gunn-Peterson 



WFC consists of a close packed mosaic of 4 thinned EEV42 
2kx4k CCDs, with a pixel size of 13.5 /um corresponding 
to 0.33" /pixel and effective field of view of 0.25deg^ The 
survey consists of single GOOsecs exposures in 5 wavebands 
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Figure 1. The gri and riz two-colour diagrams. Data from 8 
pointings (~2deg^) of the INT WFS is show. All objects are clas- 
sified as stellar and have K21.5 The newly discovered quasars 
(over plotted as squares) sit apart from the stellar locus. A model 
colour evolution track for quasars is overlaid. Filled circles indi- 
cate a 0.2 step in redshift from the larger circle at a redshift of 
4.8 in this model. 



(ugnz) over an area ~100deg^ to a nominal 5a limiting mag- 
nitude(Vega) of 23, 25, 24, 23 and 22 respectively. 

The CCD mosaic data from all survey runs i s pipeline 
processed in Cambridge (Irwin and Lewis (2000)) to pro- 



duce photometrically and astrometrically calibrated images 
(rms(internal)=0.1"), and morphologically classified merged 
multicolour catalogues. A preliminary photometric calibra- 
tion has been applied to the data to the ±0.1 magnitude 
level. For the current work we have only used data with see- 
ing of <1.67 arcsec (FWHM of 5 pixels) and stellar ellipticity 
(due to trailing) of better than <0.2 since this defines a prac- 
tical upper limit for reliable image morphological classifica- 
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Figure 2. Quasars are readily identified, even at low signal-to- 
noise, by the characteristic continuum break across the prominent 
Lyman-a emission line. The identification of low level emission 
lines due to OVI, Si IV and CIV at the correct wavelengths rela- 
tive to Lyman-a support the redshift estimation for these objects. 



tion. The current quasar candidate sample has a magnitude 
limit of i or z<~21, the practical limit for acceptable signal- 
to-noise spectra with a 4m class telescope in 1800 seconds 
of observation. 

At the time of the spectroscopic follow-up observations 
in June 2000, object catalogues for data taken between Au- 
gust 1998 and October 1999 were available. This data is 
taken from two regions of the survey. The first regio n cove rs 
the ISO ELAIS Nl field at J1610+5430 (Oliver et al. (|200C|)), 



© 0000 RAS, MNRAS 000, 000-000 



INT WFS high redshift, z>4.4, quasars. 3 



Table 1. Photometry derived from the INT WFS. 

Vega system magnitudes 

Object z RA (J2000) Dec u g r i z Mb Mab(i450) 

wfsj 224524. 2+002414 5.17 22 45 24.28 +00 24 14.6 >23.01 >24.89 23.96 21.43 21.19 -24.9 -24.3 

wfsj 161253. 1+525543 4.95 16 12 53.11 +52 55 43.5 >23.00 >25.12 22.21 20.15 19.89 -26.1 -25.6 

wfsj 224531. 0+003358 4.45 22 45 31.00 +00 33 58.6 >23.09 >24.89 21.60 20.87 20.36 -25.5 -25.0 



the second region is in the vicinity of at J2240+0000. A total 
area of 12.5 deg^ has been used. 



3 COLOUR SELECTION OF CANDIDATES 

Figure ^ shows colour-colour plots for INT WFS data in the 
griz bands with the predicted quasar tracks as a function 
of redshift overplotted. These colour tracks were computed 
assuming an underlying quasar spectrum based on a power 
law with spectral index {Si, oc z/""'^) and with an emission 
line spectru m ba sed on the composite spectrum of Vanden 
Berk et al. (2001). The absorption model for the interven- 
ing Lyman-a forest is taken from Madau(1995). The stellar 
main sequence is clearly visible as the heavily populated 
strip in the centre of the plots. 

Candidate high redshift objects are selected from the 
multi colour data by selecting all objects that obey a set of 
colour criteria. To establish these criteria the stellar locus 
in the colour-colour diagrams has been approximated using 
piecewise linear fits in the colour space. These fits are then 
used to establish the boundaries indicated in Figure |l[ 

Once selected by one of these colour criteria all 5 ugriz 
images of an object are then visually examined to check the 
validity of the photometry. This step is required to identify 
spurious sources such as objects in the difi'raction spikes of 
bright stars or satellites and asteroids that have moved sig- 
nificantly between observations. Candidate z>5 quasars are 
selected using the riz diagram whereas 4<z<5 candidates 
are selected from the gri colour diagram. At the current 
time the foUowup of the lower redshift sample is incomplete 
due to limited spectroscopic time. 



4 SPECTROSCOPIC OBSERVATIONS 

Spectroscopic observations were obtained in two nights 
(29th and 30th of June 2000) using the red arm of the ISIS 
dual beam spectrograph on the 4.2 meter William Hershel 
Telescope on La Palma. The R158R grating was used with 
the TeK4 CCD giving a spectral resolution of 2.90A per 
pixel and wavelength coverage of 2970A. Initially a central 
wavelength of 7900A was selected but during observations it 
became apparent that a useful spectrum could not be taken 
of faint objects longward of 9000A due to strong night sky 
line emission and so a central wavelength of 7500A was sub- 
sequently used. 

The slit width was matched to the seeing which varied 
over the course of the observing run in the range 0.8"-1.2". 
During the periods of poorer seeing(>l") brighter candi- 
date objects were observed. These generally proved to be 
low mass late M and early L stars as expected. Condi- 
tions were photometric throughout the run and observations 



of the spectrophotometric standards HZ 44, Ross 640 and 
LDS 749B (ING identifiers spl321+363, spl626+368 and 
sp2129+000) were used to fiux calibrate the data. Wave- 
length calibration was carried out using a mixture of arc 
spectra and night sky lines. 

Each observation was taken as a series of 900 second ex- 
posures and coadded to improve cosmic ray rejection. Total 
exposure times range from 900-2700 seconds. Several candi- 
date objects could be identified as late M and early L dwarf 
stars after only a single 900 second exposure. Observation 
of these objects was truncated reducing the time spent ob- 
serving non-quasar candidates. 



5 CONFIRMED HIGH REDSHIFT QUASARS 

The identification spectra of the confirmed high redshift 
(z>4.4) quasars are shown in Figure ^ and the emission line 
properties are presented in Table ^. Absolute B(Vega) band 
magnitudes and AB magnitudes at a rest-frame wavelength 
of 1450A have been determined from the z band magnitudes 
assuming a power law spectral index of -0.5. We use the z 
band since it is unaffected by the Lyman-a forest or any 
strong emission lines. 

Redshift determination for high redshift objects is com- 
plicated since in many spectra the only strong line visible 
is Lyman-a(1215.7A). This is often blended on the red side 
with Nv emission(1240.lA) in quasar spectra and has an 
asymmetric profile due to absorption of the blue wing of the 
line by the Lyman-a forest. The onset of this absorption of 
the blue wing of the line is quite abrupt and we take this as 
our redshift estimate for the quasars. 

wfsj224524. 2+002414 At z=5.17, this is currently 
the fourth highest redshift quasar known. This faint ob- 
ject is identified by the strong Lyman-a line at 7498A and 
the break in the continuum across the line due to Lyman-a 
forest absorption. Lyman-a shows a highly irregular profile 
with strong absorption features within the line. A strong 
absorption feature is also at the expected wavelength of N V 
emission redward of Lyman-a. 

wfsjl61253. 1+525543 This object shows a broad 
emission feature interpreted as a blend of Lyman-a and N V. 
The wavelength of any associated C IV emission line is in a 
region of the spectrum heavily effected by emission from the 
night sky making detection of an emission feature difficult. 
At z=4.95 this quasar is the seventh highest redshift quasar 
known. 

wfsj224531. 0+003358 Strong Lyman-a emission 
dominates the spectra of this object. Two clear absorption 
features, characteristic of the class of Broad Absorption line 
(BAL) quasars, are seen associated with N V and C IV emis- 
sion lines. 
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Table 2. Emission line properties 

Object z Lyman-« C IV 

Peak (z) Cent.(z) EW obs EWrest Cent.{z) EW obs EWrest 

wfsj2245+0024 5.17 7498A (5.17) 7450A (5.16) 594A 115A 

wfsjl612+5255 4.95 7237A (4.95) 7318A (5.00) 470A 79A 9112A (4.88) 308A 52A 

wfsj2245+0033 4.45 6628A (4.45) 6634A (4.46) 276A 5lA 8400A (4.42) 
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Figure 3. The histogram of Lyman-a equivalent width for SDSS 
quasars with redshift z>4.4 is shown with the measurements from 
the three z>4.4 INT WFS quasars over plotted for comparison. 



The FIRST and NVSS radio surveys (Becker et 
al. ( |l995| ); Condon et al. ( |l998| )) cover our survey regions 
and no radio emission with S20cm>lniJy(5cr) is detected for 



these quasars. There is also no evidence for X-ray emission 
from ROSAT All Sky Survey or PSPC pointed phase obser- 
vations. 

Figure ^ shows the equivalent width distribution of 
Lyman-Q for the SDSS z>4.4 quasars compared with our 
new quasars. Despite the strength of the line in our z=5.17 
quasar it is comparable with the strongest lined objects in 
the SDSS sample. 



6 DISCUSSION AND CONCLUSIONS 

The majori ty of the ^^250 z>4 quasars that ar e currently 
known (see ^ttp : / / www. ast .cam. ac . uk/ ~ quasars ) have been 
discovered using photographic plates but despite the rela- 
tive ease with which z >4 qu asars can be discovered (e.g. 
Storrie-Lombardi et al. ( |2000 ) and refs therein) for almost a 
decade the highest redshift object known was the optically 
selected quasar PC1247-I-3406 with z=4.89 (Schneider etal, 
1991). The addition of the z band to multi-colour surveys 
has changed this. The fi rst quasars with z>4.9 was discov- 
ered by Fan et al. ( |l99c| ) using SDSS data in the griz bands. 
In this paper we report two new quasars at z>4.9 bringing 
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Figure 4. An estimate of the survey completeness is derived 
based on the model quasar spectrum described in the text. The 
fractional completeness is shown as a function of redshift The con- 
firmed quasars are marked as black circles. At the current limiting 
magnitude the survey is not sensitive to photometric errors and 
so there is little change in completeness with magnitude. 



the total sample of published quasars with redshift z>4.9 to 
nine, of which 6 have z>5 (see Table 

We use our sample to make the first estimate of the 
space density of low luminosity quasars at z~5. Figure ^ 
shows the completeness of our survey with redshift. At the 
current magnitude limit for the spectroscopic sample the 
survey is not sensitive to photometric errors and so there is 
no magnitude effect on completeness. This has been com- 
puted using the approach described in Fan et al. (2001) 
where we use the Madau(1995) IGM model and an intrin- 
sic unabsorbed quasar spectrum with a power law contin- 
uum spectrum with a gaussian distribution of spectral in- 
dices (mean=— 0.5; dispersion(a=0.3) and an contribution 
from emission lines. The relative emission lines strengths are 
taken from the S DSS composite quasar spectrum (Vanden 
Berk et al. (2001)) normalised to a Lyman-a line strength 
distribution (Rest frame Equivalent Width: Gaussian with 
mean 92.91 and sigma 24.0) from Wilkes et al. (1986). Using 
this completeness we calculate the effective volume assuming 
Ho =50 km s"^ Mpc-\ qo=0.5. 

The INT WFS is ongoing and to date ~10% of the fi- 
nal survey imaging data has been investigated for candidate 
quasars. The final survey is expected to contain 10-20 high 
redshift (z>5) low luminosity quasars (i<21.5 and should 
therefore constrain the luminosity density to ~20%. This 
statistically significant sample of low luminosity quasars will 
allow the determination of the faint end slope of the quasar 
luminosity function and the break point in the distribution 
(analogue to L* for galaxies). This break point is required to 
allow estimation of the contribution of quasars to the ioniz- 
ing UV background at high redshift. Our results, when com- 
bined with shallower surveys such as the SDSS, will allow 
detailed studies of the AGN phenomenon at high redshift 
over a range of UV luminosities. Our knowledge of the rela- 
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Table 3. Known High redshift (z>4.9) quasars 



Name 



Mag 



Refs: 1. This work : 2. Fan et al. (20001:); 3. Stern et ahJ^gQCl); 4, 



Zheng et al. (2000); 5. Fan et al. {2000a); 6. Fan et al. (1999); 



tionship between accretion and the formation of black holes 
and galaxies will then not rely purely on the properties of 
the rare extreme luminosity quasars. 
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Figure 5. Finding charts are taken from i band images. The 
charts are 4.4 arcminutes (800 pixels) square with the quasar 
marked by the cross hair. 
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